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ABSTRACT 


A watronal method 1s presented, to predict the penetration 
of free-falling objects anto’ deep-sea sediments by, combining 
peoven) cmpiricaltheoenives: tromrthe field of sell mechanies: Wath 
known hydrodynamic phenomena. The impact velocity of the ob- 
mecte and the shear strength profile and density of the sedi— 
Ment are assumed to be»known.: The penetration problem was 
porved through the use jot (a computer by equating Ehe work 
meme Ciring penetration to therenergy of the object falling 
through air and impacting onto a modeled: deep-sea sediment. 
THe objects were simple geometric shapes ranging in weight 
room 200 to more than 1,000 pounds. The impact velocities 
Pomeod from Zero tO twenty fLeet=-per=-second. The cesults are 
Compared -with full scale tests and recommendations are made 
Pome meeond tye method tO a walter—Sediment intermtace. pine 
metnod successfully predicts the penetration o£ objects into 
weak, Saturated, sediments within the accuracy of the state- 
Of—the-art techniques for measuring the sediment. mechanical 
properties. The impact duration time was observed to be 
Betakively) constant fandlindepencent of object. velocity,. shape; 
Priaeweloght implying that at:may be a unique property of the 


Gyramic behavior of a sediment type. 
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I. INTRODUCTION 


A. PURPOSE OF INVESTIGATION 

Various types of objects will be Lost whieh fall senrougn 
the water column and impact on the ocean bottom aS man con- 
Eines to move more freely on the Surface and through the 
oceans. Such objects can range from small and very expen- 
Mice nSLiument packages Co Eleet Submarines. Lt Ls important 
tO know exactly how far they have penetrated into the bottom 
Siermg Search operations Since complete penetration below the 
Pwmerace’ Of the Sediments may require a search by other than 
wieweal Or acoustic methods. 

gmie welocity andjfinree fall attitude of an object that has 
been located and photographed on the bottom may be estimated 
meoMmcie penetration distance if the physical .and mechanical 
Properties of the sediment are known. The impact, velocity 
Pubmacettuce could be potentially usetul, in meconstruecring 
the event and also for determining the amount of energy trans- 
ferred to the sediment and radiated from the point of impact. 

The recovery method and force necessary to break a lost 
object free of the bottom during a salvage operation will 
depend to some degree upon how far the object has penetrated 
ineo the bottom sediments. 

The fields of soil mechanics and foundation engineering 
have been concerned mainly with static situations including 
mie prediction of the very slow settlement he a. Strueture 


undergoes as it is erected. The weight is thus applied to 
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the soil in progressively greater amounts. Deep ocean 
bottom installations may be preassembled, floated to their 
Site, and allowed to. free-fall to a postion on thes botton-. 
The free-fall method of implantment is potentially useful 
for ocean structures since it decouples the structure from 
Ene influence of surface motion and tethering lines: )  Ppre— 
eernton Of an accurate anitial penetration depth upen ‘econ— 
Mere tiecretore, becomes a Critical ‘design factor: 

An indication Of the sediment's mechanical and physical 
Droperties canbe ebtained by observing probes of known 
Moca COntrgquration as they penctrate-the -boctom since 


Bemetration depth depends on Sediment Strength: 


Pee ccOory OF INVESTIGATION 

The purpose of this investigation was to develop a method 
Pemprcarcte the penetration depth of objects dropped antora 
Simulated deep-sea sediment. Experimental data was available 
to test and verify the analytical method that was developed. 
The full scale experiment employed large, heavy objects and 
it was necessary to let the objects fall through air to 
allow them to reach a velocity at impact equal to the highest 
expected terminal velocity of free-fall through sea water for 
very dense objects. In addition, dropping the objects in 
air eliminated free fall stability problems as the impact 
attitude was effectively controlled and the effects of 
angular momentum on the impact were minimized. The sediment 
used in the experiment was of re-worked continental origin 
characterized by low shear strength and high water content 


typical of deep-sea and pelagic deposits. 


9 


=a 





Thea 
in 





The initial boundaries and limitations that were imposed 
on the study were designed to simplify the complexities of 
Paco roplen: to, perma. whe en eee be undewstood; and 
bowcetain the ultimate objective jof. pracklealzaeppo lication 
the method at sea. The conditions are summarized below. 

i... The impact,velocity at the. top. of, Ehensedimene,comunn 
Too the range. from-Zero €O,.20 feet—per—Ssecond, | Zero 
Peo vocity at. the top. of. the sediment column anvolved weleéacing 
the object at the sediment surface and its acceleration 
mamough the soft woper.layers before. coming to. rest. ~The 
Heper limit of 20 feet per second was believed to be a real- 
Pomc ecerminal-velocity, in sea water £on most, objects: 

2. Free fall was assumed, or gravity was the only down= 
Werder orce acting on the object. 

Palas OMe eCuEonapesewere Limited Lo simple eeoneeeuce] 
Weriecotablashed hydrodynamic properties from theoretical 
Sateculations or empirical tests:, Configurations designed to 
Beneerave, such aS.a,.coring tool. with a dong barrel end barge 
mass On one end were not considered. 

4. The sediments were taken to be cohesive, with no 
engle Of internal friction. . They were also assumed to tbe 
low strength, less than four pounds per square inch, and to 
have high water. contents typical of marine sediments. 

be, The Sbhudy. di d..nok, mc hideythe Sftecks, of janyieseeile— 
Pemeaatteh Lottie pencirauion.. 

6. The minimum factors necessary to predict penetration 
dépth were considered to be: the impact velocity; geometry 
and mass of the object; and the shear strength profile and 


Gepctay On tne Sedament . 10 
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G2 HISTORICAL BACKGROUND 

The loss of the U.S.S. THRESHER and the problems involved 
im the Location of the hull which had penetrated deeply inte 
very soft) sediments served; to emphasize: the fact that diceve 
was known about physical phenomena associated with deep-sea 
sediments and the mechanism of penetration. This disaster 
Prompted the original full scale experiments carried: ouikeby 
tie Naval Civil Engineering Laboratory and used toe verify 
the method developed in this study. The THRESHER disaster 
elise pointed owvt—the difficulties involved in sampling sedi- 
ments at deep ocean water depths and using only a few samples 
Peoeeqd ina Laboratory <o predict as precisely as) possibile 
the sediment's dynamic behavior in-situ. This study is only 
eeseep toward understanding all the variables.and their 
fetaeronship to each other. 

The following is a review of the state-of—knowledge before 
Enevtnvestigation and the technical background. The variables 
involved and the development of other penetration equations 
eeerartscussed . 

i. Stave-ot—Knowledge 

The problem of earth penetration by projectiles is a 
classical one of terminal ballistics and an immense amount of 
Dotw theoretical and experimental literature exists on the 
Subject. Empirical attempts to determine the thickness of 
earth embankments needed to protect soldiers were made as 
panly as L742 TYoung - L969): | Im recent years: ihe sree-fae 


penetration problem has been studied to provide protection 


ei 








against artillery, to estimate forces acting yon ayseructune 
during earth impact such as would be required by orbiting 
nuclear power plants that would have to remain intact during 
re-entry and landing on earth, and to determine physical soil 
BeGperties remotely:-by an acceleration-time record transmmected 
heen 2 projectile’ penetrometer. The studies have .involved 
Becestclds.of engineering, mechanics, fluad mechanies, geology, 
Soil mechanics, aerodynamics, wave propagation, as well as 
Stger disciplines. They were concerned with the problem o£ 
aepecdy talling with—-high velocity through, air and impacting 
@aenacd continental soils. Bibliographies .of -~prior investi- 
G@ebtons are available [Schmid, 1969 and Young, 1969]: Dynamic 
Seuaees in the field of soil mechanics on and -are.Sstill in 
Eeetinlntiancy., Whether they can be extended.to .structunes 

and the penetration problems in the deep oceans probably will 
Bemaiaed CucCStion for some time into the: future for,£ouxr basic 
reasons. 

First, bae velocity of “an object free-falling ehrough 
water will be much less than the velocity of a similar object 
Pamddbe | pAt very» Low velocities, £Lrom zere to 10, feet per 
Second, the inertia of the object may not be important com- 
pared to other factors. In water the added mass of entrained 
fluid must be accounted for, but probably would be.ignored in 
ekt 3 

Second, continental soils are quite different from 
deep-sea sediments. Sediments in the ocean basins typically 


have very small grain sizes, often clay sizes of .005 mm and 








smaller. The rate of sedimentation is measured fractions of 
mueLlinerers per year with very old layered materially found 
Beltacively close to the top of the Sediment column. =, peer— 

sea sediments are completely saturated with interstititial 
pore water often accounting for halt their bulk volumes “the 
More Water iS not the pure water commonly found on Vand," bue 
the chemically complicated sea water which causes flocculation 
Ste some Sediment particles and the presence of authigenic 
Himerals near the sea—-sediment. interface.” -<Lower in” the’ sedi= 
Heme column, GdlLagensis and cementation between particles 

Seeuc forming a Loose fabric often with the porosity berng 
Esbarrvyely Constant with depth. The effects of the Low 2 
degree centigrade temperature and high 6,000 pounds-—per-— 
‘mese-ineh pressure encountered at the mediam ocean depen 
create an entirely different environment than found on the 
Continents. The degree of reworking by benthic organisms 
Mayeatso be important in Getermining the mechanical “character— 
Paces Of the upper layers of the sediments. 

Third, when measuring the mechanical properties of 
the sediments, the amount of disturbance that the semple 
undergoes as a result of the sampling process, handling and 
preparation for testing, temperature and pressure changes, 
Pigmecganic Growth or GecompoSition is difficult te precisely 
estimate and may vary from one sampling and testing technique 
Bemanother [Richards, 1961]. 

Fourth, if the sediments in the upper layers, are 


treated as a suspension of particles in a fluid, which may 


eS 


sholo yley 


:, 








be a valid approach to sediment dynamics if the water con- 
tent is over 100 percent, then another problem arises. The 
fluid dynamics of non-newtonian materials have been researched 
for conduit, flows such, as» slurries in pipeldines,ajbulb) veryieeew 
investigations of external flow about a body [Pazwash and 
Robertson, 1969], and unsteady flow at a boundary have been 
made. 
we. Prior InvestigqalLions 

There arenonly a few records of objects fchat have 
been Tost and, recovered. The instances where things have 
been recovered or photographed to document and determine the 
penetration depth together with adequate data about the 
Heenantcal properties, of thes sediment are virtually non- 
Saeeont, The same Lacks of data eCxiSts in repokts written 
Boourt the few cases where penetration was investigated. The 
Pepeces typically. Set, forth) equations) wrthvempifrical constants 
and contain no information about the controlling sediment 
parameters. 

Much laterature and empirical information is available 
Bucie pile drivang. The forcing of a long piling anto Sort 
Sediments iS, however, quite different from the free-fall 


Benertration problem. 


D, TECHNICAL BACKGROUND 

All of the variables were studied, then the development 
of the previous penetration equations was researched before 
an attempt was made to solve the penetration problem. A 


review of the known variables for which data exists was made 
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to develop a method that would be of use at the present 
tame using eee data and techniques. 
1.7 (the Vebuabilies “involved 
Penetration may be a function of more than 20° a@vtfer— 
ent variables and may be described by the following equation: 


£ eer M, BS Z, ©&, 1G Ss, ay, Vos Ps! dB co) a Xe Ss, g ‘ By 
(2) 


eneGs Tent eaaése ph Cee F bee 
where, 

are —.impact velocity 

M ~ mass of the object 

A. =, €ress, sectional area ofthe: object 
Z ~— dépth of penetration 

rr -> surface roughness of the object 
Sy, pLebject, Shape 

2 es -veud ratio of the sediment 

E = Young's modulus of the sediment 
vy = Poisson's ratio of the sediment 
P, =.density of the sediment 

US i> Gene Or PSnerrat Om “EO 7oeCeur 

@ =~ mineralogic and chemical composition of 

the sediment 

Me, WuiuSeCOSiboy Gf Ehe Sediment 
Scio Serengen Of EMe  Searmenc 
ios SpetmeoDbiaty of Ene Ssecameme 

S > Sensitivity of the sediment 

w = water content of the sediment 

7 (Gi lation of the sediment 

@ = Grain size of the’ sediment 


tS 











Sry 
| 


= post depositional changes of the sediment 


p =VWGensery ob (ene water 

oe = viscosity of the water 

hy ® > wetern deren 
a. = adhesion of the sediment to the object 


Such “a isst 1s self explanatory “and “is Mot exe Rucive. 
Some of the above properties are related to others, and the 
parameters may thus be reduced in number by a dimensional 
analysis [Liu, 1969]. Some of the properties, however, are 
met amenable to-a-mathematical analysis, for example the 
Gear Size distribution’ and the grain's meoxphology and 
SerentacioOn in’ the sediment fabric. A phenomollegical 
Beeieeac to the problem may, thererore, be réquised Maga, 
feo TO put all of the variables antova sSystenVoL 
Suueetons may be an impossible task, and for thiseredson 
Prem Eredtecrion Of penetration: into deep-sea ‘sediments may 
fomaon Somewhat Of an art and’ not susceptible to%dan exact 
Solution, for some years. 

The deep-sea sediment particles accumulated over 
Mebivons of years Occur in many types“which reflect unique 
Beeesitional conditions of the earth's history. ~The 
mechanical properties vary widely, as subsequent to 
deposition the sediments may have undergone secondary 
cementation, diagenesis, recrystallization, and a degree of 
reworking by animals in the top important layers [Smith, 
1969]. Some problems involved in learning more about deep- 


Sea sediments are set forth in the following paragraphs. 
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Thedetails of the sea=sediment,.A htewsace sand uupoer: 
sediment layers in the depths of the ocean are only grossly 
known the world over. In the most ues SCceen ,, The wert 
Atlantic, Keller [1967], estimates one core has been analyzed 
Ser sediment mechanvecal propexties for, every, 30,000 sauianre 
ieees OL OCeCan floor. Less data is available in the Pacteac. 
Siesvarialion of the ecean bottom's natuxe and the high cost 
Semovout 1,000 dollars per foot of core recovery and analysis 
toa amit any Great expansion of knowledge. 

Much InpUiSto-khe design of foundations on, lend) vs 
Bese Od past Local experience incorporated.into Local 
mueraitng Codes: = On land, tests can be conducted and equations 
Geveloped that £lt special circumstances in cértain areas et 
Senaracrably less expense than in:the océans.,. Also -a past 
tiaseny exists Lor land areas Of anterest, whereas ‘compara— 
mavely little has been known of déep-sea sediments until Veny 
feet times: Although the horizontal area variability of 
Geer -sea Sediments may be less than that of continental 
Pecmmenes, tO Cnipirically Categorize the seventy percent of 
Enewearth's surtace covered by water and anticipate a return 
in less than the number of years that have elapsed since 
Structures have been built on land would be folly. 

2. The Development of Previous Penetration Equations 

A review of the development of penetration equations 
emphasizes the importance of the various constants that act 
Bencontrol the penetration depth. The discussion (that; fel— 


lows combines the outlines given by Schmid [1969], and 
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Young [1969]. First and most basic from. Newton ts second 
law, 


Ete) Met) 2) 


ie Can be shown that penetration depth is deétermmined it 
Merocity 1s know “ae aj function of time. " Tn @eneral thas as 
Meeeche Case Lor unsteady motion’ at a boundary. 

Other eclassical.equations are also based on Newton's 
second law where resistance to penetration is a function of 
Miemunsctantancous velocity, WV; 

aM iC 


Mer ieg tapi aa Gil eee fo~ (Wa) (3) 


Ween postulates, that during the first phase of penetration, 
at high velocity, the motion is governed by resistive forces 
Peepertional to velocity to Some power, and at the final 
Stages of penetration the resistance is independent of velo- 
ewey. Lf the functions, fie fo, pe ee are takén to be 
SCemstrants and the highest power of n assumed to be 2, then 
equation (3) Simplifies to; 


=—Vx = a+ By + eR . (4) 


Mieke lat high velocities the forces controlling penetration 
aaa epLoportional to o Or analogous to drag force ain flume 
slow. At moderate velocities resistance will be proportional 
Pomene first power of V,.analeogous to viscous resistance, and 
eee last Stage, proportional to a constant foree indepen— 
Gene Of velocity. According to how the constants, in) equation 
(4) are treated, the following classical equations can be 


derived after integrating twice with the initial conditions, 
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WA ee Vi, at tite. OF Land final -~conditionsnwVas Lolo ete 


mv? 
Robins—Buler ) 2.5 a (GonStant Kréesa sting! force) (5) 
* 
M VG 
Poncelet we ag 3 in (1 * — =) (resisting forces (6) 
VY a a Rees 
constant and Similar 
to fluid drag) 
aN : 
Resal eee oe (1 + —~) (resisting forces C7) 
(64 


Simitan te £haae 
ViESCOSLEyY and 
fluid drag) 


fie sevuestion relative to the use of the above equations Kests 
Weelecne COnStants, a, 8B, y, which are determined by the gqeo— 
Metry and mass of the object and the type of medium being 
pemetrated. Schmid [1969] discussed several empirical 
Pee eeol1 On Gquatzons from the iiterature whach ace varna— 
Seuens Of the above classical equations. The constants were 
aetermined from high velocity impacts on continental soils. 
Thomason and others ~LI6S8] used Poncever ss equarmiven 
boraetermine the» amount’ of Charge required for an explosive 


anchor system. 





2 a ms 
Zz 3b iiras COW = Vo) (8) 


where 


p = W = weight of projeeubke 
A normal frontal area 


a = constant related to the shatter strength 
of pthermedium., (for + "mud" ya- 620) pei) 

b — an inertial coefficient foe the material 

woe a Lorn. factor, Usually i126 








Evaluation of a and b is a problem and unfortunately 
aecomparison of \predieted and actual penetration was not 
given. However, the penetration depth's dependence on high 
welocity was illustrated. 

Young, [1969] presented a’ semi-empirical solution te 
the penetration problem derived from 160 full scale tests on 
Bencanental soils with. impact velocities ranging from 100 to 
#7000 feet-per-second. Penetration was treated as three 
M@idses: (1) ampact, (2) moving through the soil, and (3) 
Peng to rest.._The phenomena of entry is not the same’ as 
Boemmeng, through the soil, sbut the surface effect of Entry was 
Bema. GO be negligible after the object had penetrated three 
BeavecLaneters. The basic form of Young’ S equation was: 
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where 


n 
| 


= a Shape parameter of the ebjecte 

A > Gross sectional aréa ‘Of the ‘ob ect 
Wo Weight of ‘the object 

V, ~ impact velocity 


N = constant dependent on soil properties 
averaged over penetration distance 


Two equations were determined for two ranges of velocity, 


one for less than 200 feet per second and one for greater 


than 200 feet per second. The low velocity equation was: 
- 2 
re alae Co ae a ae AS) (10) 
and nomograms were given for easy solution. All Young's data 


fitted within + 20% deviation. Again it should be noted that 
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the soil property constant, N, could vary by 2 -‘fector rom 
to 50 and “was arbitrarily assigned by those knowledgeable 
im ithe field". 

Lf» penetration aswtreated as being resisted by yuccoms 
forces only, the,followang results according to) Schmid [1969]; 


Mx > A Wide) (la) 


which has as its solution 


ve eRe 
ES Fat (Es ) Ci) 
a 
where 
Au 
ctl re 
When t goes to infinity, the maximum penetration becomes: 
nS 
Z > — (14) 
a 


but again resolving the viScous parameter of the sediment, ym, 
ema eestablishing a realistic time, t, become difficult: 

Soils do not behave as isotropic, linearly elastic 
Mideemtats Cxcept possibly at very small strains [D/Appollonia 
and Lambe, 1970]. Penetration into a material causes the 
Pacer limit of the material to be exceeded and Ltheorsressot 
Siacercrty Cannot be applied directly to the problem. 
tmeories of elasticity and plasticity applied to, problems 
GE déformation assume idealized stress strain laws which are 
not well suited to soils and more basically ignore the 
important dimension of time, and assume that the state, os 
Serain iS instantaneous. 

Various viscoelastic models, Figure 1, with combina- 
tions of Kelvin and Maxwell models, of a sediment can be 
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assembled and solved by Laplace or Fourier transforms, if 


time dependence is linear. 


elastic rebound 


_ permanent deformation no permanent 
deformation deformation 


ae 


MAXWELL KELVIN 


Figure 1]. KELVIN AND MAXWELL VISCOELASTIC MODELS 


However the problem of determining the necessary soil param-—- 
Saeno 15 NOn-trivial, and to’ get a unique Solution; several 
Sema Of parameters are required [Schmid, 1966]. 

Modern advances in rheology and the continued develop- 
feaceor the finwte element method with appropriate yield 
criteria such as developed by Noh and Wilkins in 1964, and 
Saewaned by Schmid [1969], and the analysis of the anitial 
settlement problem by D'Appolonia and Lambe, [1970], may be 
Weesut but again the problem of relating all parameters and 
determining the stress-strain behavior and yield criteria of 
a Ehneologically complicated material such as deep-sea sedi- 
ments, may remain elusive for the near future. 

In very weak, high water content, deep-ocean sedi- 
ments, the penetration problem is somewhat similar to the 
problem of ship Slamming. Studies of ship slamming and 
unsteady hydrodynamic impacts at a water surface have mainly 
been concerned with peak pressure, velocity-time relation— 


ships as a function of the form of a structure as it penetrates 








a fluid boundary and not penetration depth. Extensive 
Pevenature and modded Gest.records .cxist. Ocha. Lo6oi 

A methed.to estimate the impact velocity srensehe 
observed penetration depth was developed by Mandl and Givens 
[1964] after an aixcraft was lost and crashed into soft clays 
goenorthern Canada... The method developed and explicity 
solved an equation of motion for a sphere and cylinder 
Penetrating at high. velocities into a Leda clay 


Sn) EM) sess Mg a= eP=D (15) 


where 


r=" ther static strength tof “the Clay, a fumcr1on 
Of Penetration depth only 


D = a drag force considered a. function of 
penetratvon depth, velocity and” tame 


Experimental tests of high velocity impacts into clay, 
Glas than 800 feet per second) confirmed that the method 
overpredicted penetration depth. The over-prediction was 
expected Since viscous dissipation and elastic wave initi- 
eatLon were ignored. 

BEEGE a detalled examination of mathematical methods 
avawlable for analysis of unsteady motions” of an object 
Grossing an air—-water interface, Moran [1965], concluded that 
meme ot them could yield a uniformly valid appreximation to 
Eade sOlutLiOn Of surface crossing and that the best hope for 
obtaining a reliable estimate was with numerical analysis. 

Ss. Known Varrables for Which’ Data Exists 
Since all of the methods reviewed rely totally on the 


ability of the selected variables and constants to correctly 
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parameterize the penetration problem, at was necessary to 
review what is known about the variables. The physical 
properties of the object are assumed to be known but the 
Geseription of the sediment strength is ‘subject to the 
interpretation of empirical tests or, in some instances, 
phenomenological observations. Even though deep-sea sedi- 
Memes are very ditferent from continental soiis; a review" GE 
what has been learned about the mechanical behavior of fine- 
grained clays on land and what is known about deep-sea sedi- 
Mentes 1S necessary to assemble a workable solution.” Also 
the problems of sampling and testing at sea cannot be 
ignored. 

infLOrmacron about the following parameters, dence, 
Mater Content, void ratio, Shear strength, and. Sensitivity, 
which describe the ben | Sepaennts oy of deep-sea sediments, 
Peer CxiLStL in the Literature or may be easily determmned 
through laboratory tests which have been standarized and 
soon may be routinely performed at sea on oceanographic 
Wessels. 

The bulk density of deep-sea sediments is determined 
from routine measurement of the weight of a known volume of 
sediment. Typical values range from 74 to 125 pounds per 
cubic foot with the most common sediments in the range from 
feeco 109 pounds per Cubic foot [Keller, 1967]. ~The bulk 
Gensities of North Pacific sediments are generally lower 
Eman those of the North Atlantic. 

The void ratio and porosity are a measure of the 


volume of sediment in a sample compared to the volume of 


ig ae 
jionar Pr 
)  Abcoly i 


er: a 
ee titi 





| 
' 





voids, which are filled with sea water in the case of deep- 
Sea. sediments. . They are also a measure Of paertrela pone co 
point contacts and bonds [Schmid, 1969]. The measurement is 
easily made with a pycnometer that measures that volume of a 
-nown weight of oven-dried sediment. 

SCnSPELVECY, SENS Lakio' OL Ehe original pundsSseurbed 
Shear Strength to a réemolded and disturbed shear strength, 
Zmimportant (to the penetration problem sinee “thie “object 
Memes cChLough -dusturbed Sediment. during penetration. /-ites 
Piegmtance 15 hacd to estimate Sinee most Samples are idve— 
Eurbed to some’ degqree. Mitchell and Houston [1969] indicate 
Eeaeeoin the mange Of ‘sensitivity from ae eS) “typical tor 
deep-sea sediments, the mechanism that causes sensitivity is 
chemical cementation, which is not well understood by out 
mechenists but obviously pera ce 

For a cohesive sediment, the strength iS a measure Of 
Ene sediment's ability to resist mechanical deformation. 
ieage: Masses of cohesive soils generally fail along aysur- 
face and hence the strength of a cohesive soil is defined as 
Syear Strength. For purely cohesive sediments undergoing 
large deformations, particle to particle interactions which 
Homldi "cause an “angle of internal friction" and increase in 
Pies strength of the sediment are generally assumed not to 
exist. However, aS more sophisticated testing techniques 
are developed, this may prove to be an erroneous assumption. 

Typical values of shear strengths were given by Keller 


hl967)) to range: from 25 ‘to’ 2¥5 pounds-per—square-ine,, witha 
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the most common range between .5 and 1.5 pounds-per-square- 
aimch. .Thée.North Atlantic basin sediments have silliqhnely 
bagher strengths than the North, Pacific. The presence tor 
carbonates will increase shear strength, but few carbonates 
ares round below. 12,000 feet. In -thesGulf® of Mexico, Bryant 
mmeeochers [1967], found shear"strengths of: silty Clay 
bueordi tes ~Ot.the abyssal plains to range from,.3 to 1.1 
pounds-per-square-inch with the pelagic sediments of the 
Continental slopes having strengths up to 5.9 pounds-per- 
Sate re—inch . 

Richards [1967] predicted that sasxeduction of in-Scieu 
strength would be caused by sampling and laboratory testing 
Pamuiaat the reduction might be proportional, to the, grain 
suze, the amount of gases in, solution, and temperature in- 
creases that cause increased Baerounal aceELon. However, 
Paletes LeSted in-situ at S00 meters with a-vene device 
Seeeelaced well with eests trom corésitaken from Ehe Same 
area a year earlier and tested with a fall cone and laboratory 
wamemikichards, 1967]... Recent experience indicates this may 
Pevemboen an 1SOlated instance of agreement. 

Laboratory tests to determine the effects of tempera— 
ture and pressure on shear strength and consolidation have 
been carried out with greatly disturbed, remolded samples 
[Vey and Nelson, 1966]. The results were somewhat incon- 
clusive due to the small number of samples and difficulties 
in instrumentation. The results indicated that there was a 


decrease in direct shear strength with increased environmental 
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pressure for a Line grained sediment with hioh vor eaeter 
The vane shear tests with increased environmental pressure 


snowed an jpmcrease in Shear strength for plastiersediments 
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tests of (for formal sen [1970]. 









2127 
) . 
2 Dept, Of Strain 
J Note: 
ee, bad 1] esas Processing rouE ontributes 
/ u 
to the ¥ excess peace Cards & ie 
219 : ie Of thesi. i DE VeErY weak, 
thes; ET One eG 
saturat esis, together ae ae Ene Upper 
is 


2124 






few fe measured by 





Orts s : 
a Vari BIR Processtnogetion (for l shear, 
Note; 
uncon Pape Te Manual Letter 135 nd cone 
. F) 


penet jes and some 





are r ts than others. 


so ayers Or 





sedi problem are 
Or cep ‘S and must often 
be physically contained durittig~. Ling: This pré- 


eludes the use of an unconfined compression Yeeic Or any test 
Pate Kequires trimming the sample before testing. The vane 
snear and fall cone tests can be done on a Sample in its 
Original container and the vane shear has the added advan- 
tage of easy adaptation to in-situ testing. Kravitz [1970] 
compared the results of Wykham-Farrance vane, Torvane, and 


fall cone tests on weak sediments and found that the results 
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from the fall cone were more repeatable than the Wykham- 
Farrance-vane or the Torvane which Showed least repeatability. 
These results may have been expected since an important 

Bae Lor in. vane shear testing 1s 0°thel rate of beading vancdmunc 
=epardtLoOn OL°C Static and) dynamic resistances [rouse a0. 

The vane) test’ svadvantages" are= | UFiest, (LE leantbe 
Conducted in-situ and a strength-depth profile is easily 
Seacineaawheresall other methods require reentry (of ayhoLe, 
mad may be vernyvdifficult-in the~deepyecean, or teouserere 
fempling which disturbs the, sediments. “Second, Leras ta 
Simple inexpensive technique when compared to sampling and 
faperatcory) testing. “Third, it does not.-require al System ef 
Meaggnts Or free-falling objects to cause the soil tovftail 
and therefore has the potential to be used at sea on a 
Polling, pitching Ship. 

The vane test limitations can be summarized as: 
Riese, tt Can only be used in uniform, saturated, cohesive 
sediment. Second, when done in-Situ it does not obtain a 
Pants forlevaluation. Third, it imposes -a failuce suriace 
which may not be the weakest surface or relevant to the 
poeoblem being studied. 

Eden [1966] “found that the vane*wrll yieldwhagher 
Strengths with more consistent results in sensitive clay 
with shear strengths less than ten pounds per square inch 
Prem compared to a fall cone. This is in conflict with Ene 
more recent tests of Kravitz and emphasizes the problems in 


testing. Kravitz made statistically significant tests ona 
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laboratory prepared montmorillonite sample and Eden tested 
natural clays: 

Shear strength values from a vane were found to be 
15% higher than those obtained by an unconfined compression 
test on Gulf of Mexico slope sediments [Morelock, 1969]. 
tis may be partially balanced by the findings of Inderbitzen 
[1970] who found that laboratory vane measurements were 22% 
lower than in-situ meaSurements made from DEEP QUEST at a 
Sepeh soft 3,500 feet. This is in rough agreement with the 
estimation made by Crisp [1968], of a decrease in shear 
Becength due sto. an expansion of pore water ‘caused by tempera— 
ture and pressure changes. Crisp estimated that an expansion 
Siepore water by two percent would occur by raising a sample 
from; i2;,000 feet deep and 2 degree centigrade temperature to 
Surface temperature and Sneeeiiee This expansion could 
Sienrricantily alter the strength properties of a sediment 
aneeerasp Showed’ shear strength should decrease, 10et0 US 
PeTeent: 

The in-Situ tests conducted by Inderbitzen showed a 
uniform increase in strength with depth in the sediment, but 
the laboratory tests were erratic, suggesting tha ee dit ier 
ence may have been caused by disturbance during sampling and 
preparation for testing. 

Penetrometers offer the great advantage of potentially 
being self-contained instruments. They either record or 
transmit a force versus time curve from which depth of 


penetration by integration and soil properties (can, be 
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inferred.» Inv,order ,bo,extrapolate tosdifferent shaseceand 
velocities, however, the curve's shape must be a character- 
istic of-the, sediment parameters.. There are real problems 
with instrumentation and directional stabilization as the 
MPeEMeLLOMeter penetrates Since Centric impact rarely Mcemes 
Wreh other than:Symmetrical objects [Schmid, 1966]. Three 
Sener problems also complicate the penetrometer approach to 
Secermining sediment properties in the oceans. First, te 

get deep information, the penetrometer must have enough 
energy to penetrate, and at low terminal velocities in water 
Paewels probably require. a smooth, known, driving foree. 
Second, viscous forces may predominate on small penetrometers 
Homma Extrapolation to larger objects very difficult. Third, 
Peneerometers do not retrieve samples for analysis of other 
sediment parameters. 

Acoustic sounding and seismic methods do not discrim-— 
foams ecnough Of the structure of the upper layers of sedi 
Memeaseimportant to the penetration problem to allow theix use 
Deerenoresent time. In the future the use of hagh frequency 
pulsed sources mounted on deep vehicles, hold promise for 
Reeeer resolution of the upper layers. 
| The methods of determining sediment strength that 
Mewes been discussed wére developed for use in classical 
foundation engineering where a safety factor based on 
engineering judgement would be applied to the final result. 
The penetration problem is quite different as it requires an 


Sxact answer since an’ error factor of two could mean the 
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difference between an object being completely buried below 
the surface or having one half of the object exposed and 


easily located. 
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iL: METHOD .OF ANALYSIS 


Prior to a discussion of the method used to solve the 
penetration problem in the air—to-sediment case and the 
modutLications required, for use at sea, it iS Suitable to 
quote directly from Terzaghi and Peck [1948]. 

"Because of the unavoidable undertainties involved 
in the fundamental assumptions of theories and in the 
numerEcal Values of the soil constants, Simplicity as 
Stomuchy Greacer ampertance than accuracy. Tf a theory 

Zroe Simple, one’ can readily judge the practical con- 

Sequence Of Various conceivable deviations. from the 

aesumpcitons and act accordingly. Lf a theory is com- 

Ppercaced, it Serves no practical use until the resulles 

aa Condensed into graphs or tables that permit rapid 

Pwerdacton Of the final Gquadtvons On the basis of 


several different assumptions." 


A. ASSUMPTIONS THAT PERTAIN TO THE SEDIMENT 

Deep-sea sediments are not well defined in the literature 
Sealine with Soil mechanics. Therefore a dist of amportant 
assumptions considered in developing the penetration pre- 
“diction method are stated: 

1. The sediment is completely saturated with sea water. 

2 The water content of the sediment. is, 75, percent, oF 
more, as is typical of marine deposits in the deep ocean 


basins. 
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3.;° Organic material that ‘would tend to %bind ther sediment 
particles: together is either absent or is present ian very 
Small quantities. 

4. The sediment will fail along its weakest plane. 

5. The only shear resistance is cohesion. 

Gen The anitial penetration takes’ place “in: sucha’ shore 
mames that the GSEeet or permeability ws nid.” inva seudy of 
pame-dependent deformation of clays Barden [1969] noted Ehat 
emencefect Of adsorbed water on individual particles was ‘that 
Zee precscscure equals the applied pressure at the £iest 
epplication of load. 

we Any lathiettcatvon or cementation that may have ‘developed 
Belong individual particle-to-particle contacts is broken at the 
SnecsOn impact: This implies a portion of the ‘sediment is 
Bemolged “at impact. If the sartienss do" not react "With one 
another by other than mechanical means, the mixture will 
behave as a viscous fluid [Peck and others, 1953, Pazwash and 
Robertson, 1969]. 

8. Shear strength varies in an irregular manner with 
meptie The variation of velocity gradients is greater in the 
deep ocean basins than over shallow areas [Ewing and Nafe, 
1963] which implies that deep ocean sediments have been 
deposited in layers. This is confirmed by visual observation 
of deep-sea core samples which often show horizontal layers 
@e ditferent colors and textures. When studying in-situ 
Strengths, Inderbitzen [1970] found that horizontal vari- 


merous were much Vess than vertical variations. He also 
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found that the percentage of variation in shear strength 
decreased with depth, which may indicate that biological and 
chemical processes are taking place in the upper layers and 
minimal consolidation. occurs. 

9. The bulk wet density of the sediment can be considered 
constant with depth in the surficial sediments. The shear 
strength of a saturated sediment depends to a large extent on 
the - Se arrangement of the particles and Scott’ [1967] 
Meade Out thaticoils can have different Structures atthe 
Pomemoensiiy) but may not have the same structure at distrer— 
Sent densities, which implies that shear strength is not a2 
mePoece Lunction Of.density. X-ray diffraction and polarized 
Mma eenectOscope Lechniques, cannot reveal the déetaiilsvot 
Peet coves Configuration or orientation in fine-grained marine 
Sediments. However, electron microphotographs of Sediments 
feem the Gulf of Mexico show them characterized by a loose, 
open arrangement of randomly oriented particles and that 
particle rearrangement resulting from gravitational consoli- 
Betaon is not apparent [Bowles, 1969). 

HOlLe The energy transmitted to the sediments in the form” 
of elastic or gravity waves from the initial impact shock is 
eg idigabie inpeompareson to the total energy dissipated during 
Bpenetration. Impacts on soils and fluid media have been 
Meudted £or a long time, but there is a very difficult instxru- 
mentation problem with detecting energy transmitted away from 
mae impact site. According to Hunter [in Schmid, 1966] the 


Hertz impact theory, which refers to impacts in the elastic 
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Range, states that the portion.of total: emergy Ehat as 
Vibrational energy iS very small as long as the impact 
velocity is less than the velocity of soumd in the medium. 
Schmid [1966] measured the distant dynamic disturbances 
caused by a small penetrometer impacting imto materials whose 
mBieOlLogiical properties were known or easily tested, such as, 
Greases, wax, and homogeneous soils. High frequency vibration 
was detected in the hard materials and it was determined by 
use of surface accelerometers that the largest amount of 
energy propagated from the point of impact was 4.3 percent 
faoene total energy available. His calculated propagation 
Beeeqg Suggested that the motions measured, at the surface 

were Stoneley waves. Ochi and Bledsoe [1962] when conducting 
model tests of ship hull forms slamming onto a water surface 
founda that forces were very large for .01 seconds -and .pressume 
Gue to impact phenomena were over within .O7 seconds, yet the 
model's velocity decrease during this time was less than 4 
Percent. Chuang [1970] in a similar investigation, did, not 
Getect any Significant acoustic pressure variations at a 
distance away from the impact point and concluded that water - 
could be considered incompressible for impact velocities from 


Beco, SO feet sper second. 


=. METHOD OF SOLUTION FOR THE AIR-TO-MUD CASE 

A rational, phenomenological approach was used to solve 
the penetration problem rather than an empirical method 
because of the poor state of knowledge that exists relative 


to actual behavior of deep-sea sediments. With engineering 
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judgment, this method may be extended to different objects 
andedifferent cases more easily than an “empirical “approden= to 
the problem. 

is (Role of Nydrodynamics 

The sediment was treated as a homogeneous, incompressible 
fluid to enable an estimate to be made of (1) the momentum 
Eramsterred from ‘the body to the sediment, and (2) the drag 
sonees exerted on the object at high velocities. If the sedi- 
mente was a true fluid it would have an infinite strain for an 
Poaed (Stress. asTheretore, the forces: required Go anviiatre 
failure in the sediment and bring the body to rest were taken 
Boamtene —-orces calculated using the ultimate bearing capacity 
formuca developed by Prandtl and Terzaghi. 

Mreatang deep-sea: sediments as Lluids bas been an 
Bpproach applied to the study at turbidity flows and slumping 
Pere Continental Slopes.) The term “spontaneous Miquetaceion | 
is often used when referring to saturated clays and implies 
Piuaa properties. Briefly, Morgenstern [1967] states that 
liquefaction occurs when the sediment is loaded to the collapse 
stress of the sediment fabric and the load is then transferred 
feowchne interstitial pore water. The increase in pore pressure 
then produces hydraulic gradients which further disrupt the 
sediment fabric causing the solids and pore water to flow 
macroscopically as a viscous fluid. Houston and Mitchell 
[1969] describe the fluid phenomena as occurring in sensitive 
Clays which contract during shear because of their open 


fabric. If the clays are saturated, the contraction tendency 
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mesults in a,large -transfer of normal stress from «cae slay 
skeleton to the pore water causing pore pressure to increase. 
In saturated deep-sea sediments this transfer of normal stress 
mesos almost instantaneous. . Continuing Ene description ot 
Houston and Mitchell [1969], as the remoulding caused by con- 
Paming Strain becomes more pronounced, the effective normal 
meeescimay fall to avery low value.» Very Jittle -steain, is 
meamrmed sto break bonds at partiche contacts and there is 
Beme build-up of pore pressure even before a significant 
mumber of bonds are broken because of slight elastic compression 
mesthe fabric and incompressibility of the pore. water. -The 
resistance of sensitive clays due to dilational energy was 
determined to be small or absent by Houston and Mitchell [1969], 
meagiWalker [1969]. Because, first, at high void ratios,.littie 
Or no tendency toward erlaeon Gee ebserved during -sheas, reven 
Biren Strains; and second, the effective stress had been so 
reduced at high strains the mobilized resistance of the clay 
fabric was small. Mandl and Givens [1964] when conducting 
high velocity penetration tests into a clay, judged that the 
compressibility phase of impact lasted until the surface 
boundary breaks and that the static yield strength of the 
Clay was exceeded by many orders of magnitude under impact 
and the clays behaved in a dynamic manner Similar to a FL wee 
Schmid [1966] when conducting low velocity free-fall impact 
experiments on a clay also observed that liquefaction occurred: 
Another reason to treat deep-sea sediments as a fluid 


was that their bulk wet density has a narrow range of variation 
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with depth in the upper layers of sediment and the assumption 
of a median density over the range of expected penetration is 
a good approximation. 
2 UlLCimare Bearing Capacity “Equation 

The bearing. capacity of a-Soil is a*measure of Ses 
Seetrrty to Ssuppokt-a load relative to its limiting equilibrium. 
mecorcing to Harr’ [1966], the first rational ‘approach t6 
Bevetoping a method to predict ultimate load carrying capacity 
Of a Soil was provided by Terzaghi who recognized that if the 
Peece dominant factors contributing to bearing capacity, the 
meas Of Ehe Soil, the effects of surcharge, and the Strength 
Seronetrer: Of the’ soil, could be assessed, a good estimate or 
Pearing Capacity would be obtained. Terzaghi chose as a model 
Beeene? failure mechanism of the soil the solution to ‘the prob- 
fem Of penetration of Paeeaiee penne: into metal vas "developed 
iy Prandtl and Reissner. In the original solution Prandtl 
was able to neglect the effect of the weight of the metal 
moved to the influence of surface loadings on the surrounding 
Material. When the material was assumed to be isotropic and 
‘homogeneous with no angle of internal friction, Prandtl's 
solution yielded, 


dg = yD, 75.146 (16) 
On the basis of experiments and Prandtl's solution, Terzaghi 
and Peck [1948] derived the following semi-empirical equation 
Bor the bearing capacity of a circular footing. 
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where iq “"—preSsure at Surface to cause: Eamiuue (esa) 
Gum CohesPonror shear (strength) ot Ene seuiimtpe.) 
Cs Wher Weg tor eehie .qoasl 


D, (> depth of footing below the swertace «of tthe tsenl (in) 


Z 
Q 
Zi 
OQ 
Z 
= 
I 


factors which depend on the angle of internal 
fricer1on, and depth. of Looting 

The LLESt “kenm represents, a (contra butien due toute 
meee Serengtiy (of “the soil, the second the effect of material 
removed by putting the footing below the surface (equivalent 
eo buoyancy in water), and the third term, the effect of an 
angle of internal friction (assumed to be zero for deep-sea 
sediments). The equation is modified slightly for square or 
oblong footings and Meyerhof [1961] gives constants for wedge 
Pmarcone shaped footings. More complete diseussions of the 
Beating Capacity equation and slight variations in the values 
Senene “Constants can’ be found in texts on soil mechanies 
Pieuzeghicand Weck, 19437 Harr, 1966; and Hough); (1969): 
Terzaghi's equation was the closest to Prandtl's (which 
requires theoretical failure surfaces) and it has been 
Progrically verified for many years. oft was therefore 
selected as a reliable relationship for estimating the static 
strength of deep-sea sediments. The values of the constants 
and the shapes used in equation (17) are given in Table I. 

Insthe® bearing capacity equation, the shear resist= 
ance mobilized on the failure surface of any surcharge 
material is neglected, hence the equation is conservative 


and relates only to events occurring below the footing or 
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penetrating object. During penetration, static resistance 
will occur in the undisturbed portion of the sediment below 
the deepest penetrating: part, “the. other surfaces. of , the 
penetrating object are subjected to remolded sediment or 
sediment already set in motion by the penetrating object. 
gherefore the bearing capacity equation constants used. in 
penetration prediction were considered to be for shallow 
Pootings for all depths of penetration. 

Be, Une Charackeraot Failure 

Before a final method of solution was selected, 
MOEION pictures of penetrations into soft sediments and still 
Photographs of soil failures were studied. In the motion 
puecrures, the sediment was observed to move laterally outward 
from the object and have somewhat of a separation point from 
moe, object. While dropping plates vertically into a simu= 
fated sediment, Erchul [1968] found that polymer coatings to 
Bence Grag on the plates: had little effect at high entrance 
Velocities and postulated that the sediment may not have been 
Mircoueact with the plates. . He also observed that at low 
melocities the lubricants did reduce adhesion to the plate 
and considerably greater penetration occurred with lubricated 
plates Gonpared to Uniubrvcated .ones. 

To Geincea uci tative indication of sehs Gitecmser 
fluid in interstitial spaces, Rowe and others [1962] dropped 
probes into finely ground (less than 40 microns) silica dust. 
At atmospheric pressure the dust displaces laterally and 
Graters were formed. In a vacuum with no air in the inter- 


Stitial spaces, resistance to penetration was 100 percent 
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greater and the dust did not move laterally and no mounds 
or craters were formed. 

THE MOELON pPLCEUGeS- Of penetrating eob secus an sceae 
Best pit revealed Little splatter and little isuriaceudise— 
Euebance at a distance from the penetration point. ~Clark 
and Robertson [1965] observed motion pictures of a buoyant 
body exiting from water and came to the conclusion that 
Puemgy lost) due to \surface effects was small: “Moran | 1965) 
Mmus@eatrrivyed at the same conclusion. 

ay HoUAat1on_of —Mots.on 

Sanee penetration “aus a dynamic phenomena involving 
Ene Unsteady motion of a body, the first attempt at evalu- 
ecing the penetration depth involved writing an equation of 
Berean tor the object as it penetrated the sediment. — bone 
feiat solution required a Nee energy approach that allowed 
ene change in velocity and penetration depth to be determined 
moerementally with a computerized numerical scheme without 
Paercequirement to integrate the equation of motion and 
establish the instantaneous inelatvonship between ftoree, mass 
ana acceleration. 

The forces acting on an object penetrating a sediment 


Eee@eedetveri zed as a Fluid with a static strength are shown in 


mrigure 2. The equation of motion can be written as follows: 
Mala > Applied Forces (8) 
Where the applied forces are: weight of the object, the sedi- 


ment static strength, the buoyancy provided by the sediment, 
Ene force to overcome the inertia of the sediment, and the 


Grag forces. 
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Figure 2, FORCES ACTING ON A PENETRATING OBJECT 


one Change of the body Ss mOEbl1on as 1b penceravec ene 


Peement: 15 the equivalent of the resistance encountered iby 
Pac body, or the energy involved in the change is equal to 
iat expended in propelling the body. Momentum iS trans= 
Perted to a Surrounding fluid by; (1) a certain volume ef 
fluid being accelerated uniformly, (2) the fluid being put 
moeo turbulent motion, (3) compression wages ey gravity 


waves or spray [Hoerner, 1958]. All of the momentum is 
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eventually consumed by viscous friction. The last-two modes 
of momentum transfer have been assumed to be negligible for 
the velocity ranges and object sizes under consideration and 
were discussed previously. 

The, treatment. of the drag, force, or the momentum pur 
mee the fluid by turbulent,motion and viseous drag, outlined 
below, follows that given by Mandl and Givens [1964]. 

The drag forces cannot be evaluated easily for unsteady 
moron through a fluid boundary. The coefficient. of drag, can 


be ‘approximated_by,; 


e es Fis 


; e (19) 


Deo 
Zomerrst term represents the contribution due. to viscous 

G@rag; the second is that due to form drag; 6 is a constant 
depending on the shape of the body and is uSually represented 
by the first term in an infinite series (Oseen's approximation). 
Por a circular disk placed broadside to a stream B= 82 while 
mer a Sphere 6 - 24 (Lamb). The term Ch core enone to the 
Ssymptcotic case aS Reynolds number approaches infinity in 
Navier-Stokes flows, it depends on geometry and not the fluid 
Gynamic properties. Inspection of experimentally obtained 

drag coefficients versus Reynolds number relationships for 
biaurtt bodies show that Ca is a Slowly varying £unctienm at 

large Reynolds numbers. For bluff bodies with sharp edges 
Mearly all the drag is form drag. The experimentally deter- 
mined values of S38 the case of bodies with. sharp edges is 
approximently 1.0. For bluff bodies without sharp edges, 


such as spheres, the form drag also depends on the position 
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of flow separation which in turn depends on the critical 





Gange o£-Reynolds numbers. The equation cf motbtonstnen 
becomes: 
eps 1s) cee 2 
t+ = = — a ~ 
(m'*M) x Mg Ps ¥ Ax 5 A ne 5 PAC, ms (20) 


Lf the viseosity of clay Ys taken to be om the order 
Me > xX to lbs - eeey a [Whitman in Mandl and Givens, 1964], 
Ene Reynolds number for the objects that were under study with 
a three foot maximum dimension and with a velocity of ten 
fee Per second, 2S On the order of 2°x 107. Admittedly a 
Benceant Value of viscosity for a rheologically ‘complicated 
Beoosecance Such as a deep-sea clay is a gross over simplifi— 
Paeten, but the order Of magnitude of the terms tMndicates that 
mieevESCOuS GQrag Contribution to the equation Of motion may 
mentgnoced {Mandl and Givens, 1964]. Also, towards the 
Zonal phase Of motion where the effect of Viscosity is 
Prable to increase rapidly with decrease in the rate of 
Strain, the influence of viscosity will be negligible Since 
ene velocities are small and the other terms in the equation 
are larger by several orders of magnitude.: Since form drag 
represents pressure forces on the body Hee are Signi lveame 
mete snigher velocities, 1t could not be neglected: 

ae naded Mass BEEECCES 

Another reason to treat the sediment as a fluid was 
to have a method available to estimate the volume and velo- 
city of sediment displaced as the object penetrated. Womack 
[1967] verified that the amount of saturated clay moved 


Guring penetration was equal to the size of hole created. 
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The energy required to move the sediment is conserved since 
ence the (sediment is put in motion by thevebject at ware 
continue to move away from the object creating a larger 
eratcer as the object is slowed by other forces) We asi very 
important though, because the drag term of ‘theequatvon orf 
motion is dependent on velocity squared and the inertial 
Sencet of the: sediment: motion will alter the ‘velloecity—time 
fastory of penetration which has significant effects depend- 
te upon the shape and impact velocity of the object. 

Added mass is caused by hydrodynamic pressure forces 
Bare are proportienal to the instantaneous values of -accell= 
eration and according to Wendel [2950] appear ve have the 
weame value inva viscous or inviscid fluid as long as’ smali 
bodies are not involved. The added masses were determined 
ieeencing tol equations’ and methods given by Saunders [1957], 
men are tabulated in Table IL. 

The assumptions pertaining to the fluid approximatrvon 
of deep-sea sediments can be summarized as follows and are 
the same as given by Moran [1965] in a study of water entry 
ome Exit’ 

Ls Lrrorational fLlow = Viscous eriLects on the pres 
Sure felt by the body are generally small and the variations 
in inviscid forces are more important than viscous drag. 

2. Incompressible flow - impact velocities are well 
below the speed of sound. 


3. Surface Grfects are negligible. 
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The problem with determining depth of penetration 
from equation (20) .is that it cannot be integrated during 
Ebe unsteady motion of an object transiting 1a fluid boundary 
meomnone fluid EQuanother. This is true for the tollowang 
reasons: (1) m', the added mass and mass of the displaced 
Seaiment)) |S apfunctvon of depth of penetratvon for Man Vebiject 
whose submerged portion changes shape as it submerges; (2) 
eae Gominant term, that of static soil strength, varies wilth 
eae. Gepth in a highly irregular manner; (3) the location, of 
the air-to-sediment boundary condition is unknown a priori; 
and finally, (4) the flow is unsteady. 

ee “lranstormation, of the Equation te Work=Bnerogy 

fie equakvon jot mMmOtLon May Casrily be resolved sine 
an equation of energy and work which is time independent and 
Peemec1Oon Of increments Of distance. The development is 


given below. 





ae: 
er oa (21) 
dt 
Maeeaplication by dx gives 
= G_ (gx) = y GX dx 
pe eee ex a ee Mae ae) (22) 


mine first integral with respect to displacement is: 


Z 
= ax 
= / psf + 
jeg ax a2 MM Ce c (23) 


The left hand term represents the work done on M during the 
interval dx and can be evaluated if BL is a known function 
Of x. The right hand term represents the corresponding 


change in kinetic energy of the object. 
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Kinetic energy 1S a scalar quantity that depends only 
On the: mass and magnitude of the velocity. Simee the velocity 
is Squared it is always positive. The total mechanical energy 
SBeethe body is -equab: tos 
Eee oer PE (24) 
Por an object penetrating the sediment and coming to rest, 
Meemenergy of the object can be accounted for by: 


as La = 
aT Pony our mee (25) 


mars sinple balance States that the total energy E. of the 


E 
Begect at the beginning Of an interval plus any energy put 
facto the object minus any energy taken out must equal the 


final energy, E Of “Ehe Object at the end Gi the interval. 


2! 
Maem the object has Tost its energy in overcoming buoyancy, 
Beecre Soll strength, inertia of the sediment, and drag, it 
Saul Stop its motion. 

ee pateicullanr, rererence Level for snerenercates ioe among 
meeded Since energy is relative. . The potential energy 
SBepends On an arbitrary selection of datum for zero potential. 
Mmes@etme kinetic energy iS arbitrarily expressed relative to 
mie Condition at zero velocity. Hence any physical measure— 
mente Over an interval describes a change in energy. it 
follows that reference to an arbitrary datum used for the 
Pere ssion of energy will cancel. 

The work-energy method approach to solving a problem 
is powerful but dangerous because an answer is always obtained 
With no indication that the answer is correct. The method 


balances work to equal energy available and if a portion of 
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ef the energy available’ is mot considéredy om eincermreeris: 
Semnouted. works stiri donewand an answer HESULES Sin “ene 
wisteady (case of penetration, using an iterative solutaon 
ermoSst Kequires a icomputer because: Errors, in) Computation wien 
not be indicated but will appear as energy gained or lost. 
Gace answer’ will not appear wrong unless a gross error was 
made. 

in analysis Of the pencstration problem only mecianr— 
Bol ecnergy Was Considered, with the sediment and object 
Beiggd a Conservalive-system in the gravitational field. 
ihermal, chemical, electrical and atomic energies were 
meg lected. 

foe Mechod (GE Solucioen Used: with a Computer 

the numerical ateration scheme used to solve the 
megetration problem was the same for each object shape. 
meverry, it divided the object into sections or Slices in 
Eniewcame plane as the sediment surface. Then the work 
Mececoary LO Submerge Geach section and the work done by 
ewan Oot the Submerged sections, the forces due to the 
eaiettia ©: agded mass and displaced mass of the Sediment, 
Miooracc Stbengcea Of the Sediment, and drag in moving an 
Beetbemental distance equal to a Section thickness were 
Pabeulated and subtracted from the kinetic energy and change 
in potential energy of the object. Penetration depth was 
determined when the total kinetic energy of the object 
Seueleag Zero. 

A description in mathematical terms follows and a 


block outline of the computer solution is shown in Figure 5. 
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Constants and necessary approximations for each object are 
Listed in Table I. The computer programs are included as 
Appendixes C through H. 

Refer to Figure 3 for an explanation of the 


Gaertingwevons that odslow- 
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Figure 3. DEFINITIONS 


= average velocity of the object as it transits 
distanee x Ecrom position 0 to posttiom i 
(Assumes a linear velocity gradient) 


Ke 
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Vi, ~iaverage velocity of the objeci asltiel teansi es 
dastance’ x £Lrom position 1 te) pesi tion 


M > mass of the object 


my —added) mass at Section 1, pilus lithe massvor 
displaced flue? atbiseciion: wi 

Ms, = addedimass at section 2, plus-Ehe mascuor 
displaced Gluid at’ section Z 

in > m5 ms 


U ~> work done 


U, “work done against buoyancy forces as object 
b : oR: Roe. 

brans2ts LErom position 1 to poszeren) 2 
Ue > Work done against Sediment Stabic strength as 


Ene WbieCrSs Eransits from posttivon a to 
posveion 2 


U. = work-done against drag “torces as Ene object 
transits irom pesitaon Ll to position 2 


AU ~ the change an energy of the body as ai treamsived 
EsoOMm pOSteuon | Eo position 2, met ane luGiias 
change jp pOrent val energy 


Ey — total energy of the object au posrEvong® 
Ey = total, enexgy of the object at positron 2 
er Change Am potential energy Of jee Jobvjeeuiowcr 


distance x 


From equation (22) 


2 d 
U = m” vod Vv (26) 
Ln 

2 2 

Uw 72 ni" Nahe Vf 2\ m0" vi (27) 
me tollows that 
a 1 Z ‘ 2 te 1 2 - + ao 
ao = 1/2 i LA2 my, Vy VA em V5 U,, (eb Ug 

(28) 
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NOTE:: lL. Am” was accelerated from (a State candle om eZee 
velocity ‘to the velocity vat the endwot the tsection 
ever which 2t& was computed. 

2: Am’ owas assumed’ tovbe zene aibverieneworree ty had 
penetrated to its maximum dimension since 
separation of the sediment frem the Gbjeck was 
apparent ac “EhES spore. 

The incremental distance, x, is fixed, and the various masses 

ere known aS is the-work necessary to move the object from 


Pestelon one.to-pesition twos. , Velocity at. position two is 


ioknown. Therefore from the work energy, relationship 





& a 2 
— +- ee = \ 
E. Ey PE - AU Iie Mv; (29) 
Combining equation (28) with (23) and solving for Vo gives 
ane SoG ae ELIAS TAREeR Ga ce 
7 + = = é 
2 (Mi Aan + me 


i) 
Pieene 1leration scheme, when v, equals zero, ehe ‘Sum lon gene 
MmeGemenis, X, 1S the penetration depth. 

tne work done against bucyancy forces 1S a Strartgie— 
forward calculation that is discussed in elementary fluid 
Seeaanees Lexts. In the Computer iteration scheme where cre 
Sections of an object are taken to be slices with vertical 


Sides it becomes; 


UL, oR a A 
Bias )|.>) 7 A = (for all sections below the one 
being submerged) 


(for the section being submerged) 


The work done against drag forces was taken to be: 


Us =, pCy A Si x (3a) 
where y = bulk wet density of the sediment 
Cy “(drag cocetticient for the approximate ebject 
oo 


Shape that is below the sediment surface 
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A ~ the maximum cross sectional area of the object 
below the sediment surface and normal to the 
velocity tvectoxu 


v — average velocity over the distance moved, 
assuming a Linear variation an veloc. 


The work done against the sediment static strength 


was calculated from the ultimate bearing capacity equation. 


OL TuCe MN Ave x . (32) 
- where Cy =a Constant depending on footing Shape 
Na = a constant depending on the type of sediment 
ane. Looting, prokile 
c > shear strength of the sediment 


When the penetrating object had a flat plane storeing 
Paco, ehe Sediment such as a cylinder impacting on end, the 
meescs Sectional area,.A, Of the flat area and the shear 
meeemociy c, of the Sediment at Ene peneceration Level tof tne 
flat area were used in the equation above: “in he case vcs 
objects with rounded or sloping sides, a: judgment had to be 
made Of a cross sectional area and representative shear 
Merengch. The area selected was the maximum cross sectional 
Sees enornal tothe velocity,vector at the surface of the 
sediment, or the maximum cross sectional area of the object 
if the area had penetrated beneath the sediment surface. 
This method of approximating a sloping or curved surface by 
a projected horizontal area appears valid for blunt objects 
[Meyerhof, 1967]. Since the shear strength characteristics 
of deep-sea sediments in the upper layers appear to change 
radically with depth in the sediment column, a shear strength 


for equation (32) was determined by selecting an average 











depth according to the zones of plastic failure implied by 
eonationse (26) sande (47 )icco<Thisewas selected as one-half of 
the penetration depth until the maximum cross sectional area 
had been reached, then it became the depth of penetration 
minus one half the penetration depth when the maximum cross 


sectional area was coincident with the surface of the sedi- 


ment as illustrated in figure (4). 


cross sectional area, A 





a + 
ee = depth determining the value of ¢ 


cross sectional orea, A 


lia sediment surface 











depth determining the value of e¢ - 


one half the depth when the maximum cross sectional 





area was coincident with the sediment surface 


Figure 4. DETERMINING THE DEPTH OF "c" FOR NON-FLAT OBJECTS 
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The method described was adequate for curved and 


conical objects and approximates a representative strength 


in the theoretical zones of plastic failure in the sediment. 


fm alternative and more refined method of selecting A and @, 


would have been to use the differences in the cross sectional 


Bmea Erom Section to section and apply the value of ¢ at Ehe 


meen (Oh Cache S6ction.« This) increased’the complexity of the 


iteration procedure, hence was attempted only when there 


were discontinuities in the penetrating object's shape such 


Bemignt occur with a truncated cone. 


The distance of-each iteration was chosen only to 


conserve computer time and approximate the expected accuracy 


Sm tae field tests. Larger and smaller iteration intervals 


were tested for effects on the final 


depth “and, as long tas 


the interval was less than one-half of the distance between 


Shear strength measurements, convergence was always apparent. 


wend linear shear strength profile, 


the interval is limited 


muy by the desired accuracy of the solution. 


The equations. for added mass 
fertlow chart of the basic program is 
Appendix A presents a set of 


illustrate the relative influence of 


arergiven in “able shige 
Given ingrigune So: 
curves which qualitatively 


each, feree an the pene- 


tration equation. They were obtained by using the sediment 


mieat Strength profile typical of the world's deep-ocean 


basins as determined by Hoag [1970] and shown in Figure 13. 
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Circular disk moving 
normal to its plane 
or rotating about a 
diameter 








Mode of motion 


8 
m= 3 pa° 





cee 





Radius b 


Moving is 4 2 
m) =k = 
broadside Uke (Sproat) 


Prolate 


4 
ellipsoid Me =ky (4 prob’) 










Planetary ellipsoid 
or oblate spheroid 
Moving parollel to 


its polar axis 1 2 


sin? e-e,/ 1-e 


e- (sine), /1-e2 
m= pro} ere 


Elliptic disk moving 
perpendiculor fo its 
plane 


Mm, =k ( pra b) where 


] 


eo 
0.5 
! Aes eosin 26 + cos*6) d@ 
1.0 10.0 
0.637 041 0.19 0.098 


| Table Il. FORMULAS FOR CALCULATING ADDED MASS 
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PHoerner, So. &., Fluid Dynamite Drag, Practical Intommataou om 


Beteaynamic Drag and Hydrodynamic Resistanes, Chap. IL, Pub- 
lished by the author, Midland Park, New Jersey, 1958. 
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faeereactice, p. 172, John Wiley and Sons, 1948. 
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Seeundéers. Harold H., Hydrodynamics in Ship Design, sv.) ia) 
Pa +21, The Society of Naval Architects and Marine Engineers, 
2257 . 


On! 
O) 





vy . 
2 
5 
1 
; ] 
¥ 
i 
" 
" 
i 
: i 
i 
7 
[ 
& 
1 
ri 
1 
i 
t 
ss, 














Set constants for object, weight 
dimensions, coefficients, sediment 
density, iseciivon Anternvall 
















Caleylate SeEtLON, PrOopereLes or 
entire object 

1. Average projected area 
2. Volume 








Caiculate added mass for each section 
and change inaadded mass 


Calculate work done by buoyancy 





Readers Weie Impact. velocitwes 
2. sediment shear strengths 


Caletulates i. Original energy 
2. ,PoOorential” energy 


Calculate shear strength profiles 
(linear variation between points) 


| Calculate work done by sediment = 


Calculate velocity, neglecting work 
done by drag 


Calculate work done by drag | 
ogee 


Recalculate velocity using average 
yelocity andadrag 


<H verocity > 


Conpute penetration 




























Figure 5, OUTLINE OF COMPUTER PROGRAM 
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Ge; MODEFICATION OF THE “METHOD FOR USHA Ty SEA 

The method was ultimately developed for general use at 
sea with as few restrictions as possible placed on Size and 
Shape of the object and on the type of sediment. However 
Ene *only data available to test and verify the method was 
for an air-to-simulated deep-sea sediment case. 

taers arectheee relatively ‘sample alterations ethat isheula 
Bemade to| the computer program before it can be wséed- in a 
weer —-to-sediment ee thenenree jadeerat tons eshowlditachrrco 
miuerease penetration depths ‘at higher velocities and their 
mea bE Correct Wis bocineréease the importance ef ran \accurate 
Becermination ofothe shear Strength profile. 

Riese, che added masscerfects will. change by the’ ratio wer 
Ene Gensity of sea water to the density of the. sediment. This 
Gan Gasily be accomplished by dropping the unit value of sea 
eavemutron the bulk wet dénsity term in the programs. 
Additionally, the added mass of the object free-falling 
through the water must be added to the total kinetic energy 
Berenae "object at impact’ and a method for calculating the 
PaseeirOf this entrained water mass ‘taken intolaccounts 

Becona, the buoyancyterm will also» changevby sthesratio 
6f the density of the sea water to the density of the 
sediment. 

taro, i theadrag tenm wihlwchange) by the ratiosotesead 
water density to sediment density. Also a drag term must be 
Bhetuded for the part of the body that projects: above the 


Sea-sediment interface. The easiest way to accomplish this 
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would be to include a drag term for the whole object and add 
a.drag term,for. the portion of .the, object submerged) in.&the 
sediment. 

The three modifications will serve te reduce the influence 
of added mass, buoyancy and drag and cause the sediment 
Beecac) Strength, term to,become; even more dominant. thanesin 
the air-to-sediment case. Keller [1964], used the standard 
Dearing-Ccapacity equation, only, to predict the initial 
Poeeleoment- depths, of a rectangle, square, and cylinder, inte 
soft bay sediments. The results were poor and he concluded 
min pcObLems.with determining an accurate strength of the 
Sediment were the reason his calculations always underesti- 
Mated the settlement depth. 

Me correctly use the method at sea, several cores should 
be taken and their physical and mechanical properties averaged. 
eyvecaging should reduce errors caused by sampling and testing 
Breturbances. The effect of sensitivity, as it applies to 
the shear strength used in computing penetration depth, is 
Mieecrtain. Possibly the remolded shear strength is a better 
malue if it is assumed that impact shock effectively remolds 
mae Sediment to an unknown extent. In sensitive sediment, a 
range of predicted penetration depths should be calculated 
mecerding to the degree of sensitivity. 

fae orientation that the object has as it first Starts to 
penetrate has a significant effect on total penetration. 
After free-fall in several thousand feet of water, the 


meject may still be tumbling or oscillating if it is unstable, 


3) 


fenmid [L966] found that a Significant portion of a, talling 


body's energy may be in its angular. momentum which is not 











included in this method. 
For these reasons when predicting penetration depths, a 


Meomputer routine should be used to solve for a range of impact 


‘ 


Merocities, shear strength profiles, and object orientations. 
‘the Boalculations can be done by manually using larger inter— 
Beis, as Was done to check the first computer solutions, but 
since energies are being equated and the method is iterative, 


extreme caution must be used since an answer, correct or not, 


always results. 
ye 
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IV. RESULTS OF EXPERIMENTS 


A. THE PROBLEM OF MODELING IN SOILS 

Some problems with testing and modeling the dynamic 
behavior of soils are discussed before reviewing the data 
evyailable to test the penetration prediction method and the 
meagets of the tests. 

mane field of soil mechanics has traditionally been con= 
Senrned with large scale construction operations. Many of 
Seeemethods used to calculate :the stability of slopes and 
meaadation designs are based on ultimate collapse criteria 
and much of what has been learned about soil behavior was 
G@erived from studies of catastrophic large scale failures. 
Bea Construction projects are seldom undertaken without 
mepensive, full-scale, field tests of critical portions of 
_mencoundations, and modern buildings are carefully instru- 
meee EO MOnitor their settlement and reaction to seismic 
disturbances. A good indication of the state-of-the-art of 
Seale mechanics is implied in the large safety factors 
Beeonmended in texts and building codes. The lack of under- 
Beamaing is, in large measure, due to the inability to 
Parameterize and model a substance as complex as a soil. 

When the problem is narrowed to the field of saturated 
deep-sea sediments, the problems of modeling are somewhat 
Feduced. A modeling approach to investigate the dynamic 
Besponse of soils to impulsive loads was PeREMe wee by Westine 


[1966] who found that a satisfactory case could be made for 
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a model with saturated clays, but granular materials with 

angles of internal friction were impossible to Simulate, 

Since dilation and pore water pressures were hard to model. 

Prolich [1967] attempted to model.the dynamic. behavior of 

Submarine sediments, but unfortunately he passed away before 
. Beet shing his)» final nesults'. 

There are several areas where modeling dynamic pene- 

tration phenomena in deep-sea sediments may prove to be 
ampossible. 


Under dynamic loads a Saturated deep-sea sediment may 







Meehave as a bingham body.” The inertial resistance of the 

' sediment to deformation coupled with the threshhold stress, 
emamacic resistance, of the sedament, may be very. hard to 
separate if the mass and therefore, inertia, of the pene- 
meacancg body is small or on the same order of magnitude as 
the maaeie resistance and inertial resistance of the sori. 
eeemalogous problem, but Simpler by one factor, iS unSteady 
mtiuid flow at low Reynolds numbers where inertial and viscous 
mOECeS are the same order of magnitude. In the dynamic 
testing St Sots this problem can be circumvented by f£ull-— 
Meare tests using objects with large masses. 

| Mcaumilar problem involves surface effects and the 
‘adhesion of deep-sea sediments to an object. The effects 

Of adhesion have been shown to be influenced by polymer 
meetings, Surface roughness, surface chemical composition, 
eeecetera, by Erchul [1969] and other investigators studying 


mele driving on land. Surface effects may control penetration 
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Significantly if the object has a large santace area an 
relation to its total mass. A means of minimizing adhesion 
mereccs would be full scale tests with. large masses. 

me a study's purpose iis to learn more about adhesion, ‘or 
mass effects, then modeling may be the best approach. However 
Pe seaDpolation of modeled test results to objects larger by an 
meager Of magnitude may be invalid. 

Finally, when attempting to model a sediment composed of 
miay particles which are very small, to reduce them further 
Menateatltempt at modeling the interparticle attractive forces 
mea pore water viscous effects, may be impossible. Also, the 
mastruments available to measure the modeled sediment strength 
Muerte nave the required accuracy or precision for a reduced 


scale. 


=m. DATA AVAILABLE 

The Naval Ordnance Test Station, Pasadena, California, 
_eeeeerecied a large pire -from sheet piling; 3705 feet in 
@eemecer, and filled it to a depth of 16 feet with a very 
_emeseaiment from Seal Beach Lagoon, California.’ The pit 
meas used to conduct tests on the holding strength of anchors, 
meeakout forces of objects imbedded in the mud, and experi- 
tts With dynamic penetration. The penetration experiments 
_—emedrried out under. the direction of Dr. R. J. Smith in 
September, 1965. The penetration experiments are described 
below. 

fen SDeseriptiom of Objects and Tests 

Figures 6 through 11 show the six different objects 


that were dropped into the test pit. All six were fabricated 
63 
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The «sphere (was made jot ta 
heavy (steel plate rand 
weighed 561 pounds in air. 

A padeye provided an attach- 
Ment (pOLNhetoOr aL Tec imSe 





Figure 6, SPHERE 

































































Dae Noxrizontal eyirnder was made from a’.55° qallon 
orl Crum. Lt weighed 540 pounds in-aie. LE was 
metLed by a two pennant bridle with one-pennant 
mekded at each Lip. 


Figure 7, HORIZONTAL CYLINDER 
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Figure 8. 


Mae cube was fabricated from 
Beeel plate and weighed 702 
Meounds in air. It was lifted 
meyea LOur pennant bridle to 
mem corner of one face. 


Figure 9. 
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The cylinder-on-end was also 
made) froma Souigal lon iond 
drum. | For ithe £irst tests 
it weighed 540 pounds in air. 
Later a) flat, Treound plake 
was welded! over ‘the’ lito Gt 
the drum mtitche? end wthat 
penetrated, making Lt smeoth 
LO uReduce Ehe amount of air 
that might have been trapped 
there. The plate increased 
the weight to sel pounds an 
ate Tie -Gisum «was: aries. jay 


a three pennant bridle attached 


EOv its upper end. 


CY LINDER-ON-END 





























CUBE 

















The rectangular para- 
Llelepi ped weighed 

1, OL4 pounds an. aan, 
It was constructed of 
steel plate and a tted 
With, a four pennant 
brvdle to four cormiers 
of a large side. 






































Figure 10. RECTANGULAR PARALLELEPIPED 

















The cone was fabricated of 
Sse plate and weighed 





Figure Il, CONE 
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with.a, steel shell. for, structural cagsdityoand theveavyaty oF 
the shell was completely filled with fresh water for weight. 
The weight could be changed if desired, however for pene- 
eration tests they were pressed, full. at.adl times., . Dimensions 
miethe. objects are, shown, in, inches. 

The objects were suspended by a crane over preselected 
Pete tons in.the mud pat; They were released by, a; quick 
release hook that was tripped manually by pulling a line from 
Ehe side: of the pit. The objects were released at different 
merits above the mud-air interface-to allow them to.reach 
Semocities of 5, 10, 15, and. 20, feet per second as the 
Object's lowest portion contacted the sediment-air interface. 
Mae zero velocity test was considered to be a zero release 
mMeight, and the object was released with its lowest portion 
just in contact with the sediment-air interface. 

Slowsmotion. pactures, 32) frames, per second,.of almost 
Meche penetrations were available. They were carefully 
Seugied in a frame-by-frame projector to double check the 
meng Gata sheets for-height o£ drop and penetration depth, 
Bie movies also gave an indication of the stability of the 
egect aS it penetrated. Remarks about each penetration are 
maven in Table X. Times. for penetration and times for 
mebound were estimated by counting the frames of the movie 
film for cach event. The movies confirmed the observations 
Bf Chuang [1970] that air trapped below a body impacting on 
meeeuid does not alter the penetration depth.,, Any air that 


WaS trapped below the body usually vented by bubbling up 
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muangside..the object after ae had -come sco mess. FOnis7. am 
one or two cases was it observed to form a jet with a 
Significant amount of entrained mud. 

The penetration depth was determined by eae 
observation of marks painted on the objects which were scales 
Givided into one quarter inch increments. The depth refer- 
eaece was the undisturbed surface of the pit. Depth obser-— 
Meaetons were made at each corner of the rectangular objects 
end at a diameter of the axisymmetrical objects. The pene- 
Beation depth was-considered to be the numerical mean of the 
memeral measurements. Any exceptions to this method and 
meetes Of tilt at the object's final rest position are ‘noted 
mmeme remarks column of Table X. The exceptions occurred 
When the cylinders-on-end tipped completely over after coming 
Seerest. The penetration before tipping was used if it 
meld be read from the movies with an accuracy of one half 
each . 

Peeoescription oer the Sediment 

There was a real problem finding a suitable material 
=O Simulate a true deep-sea sediment, compounded by the 700 
_—_eovacrds of material required to fill the test pit. The 
Material selected was lagoonal mud, described as "“unconsoli- 
mated grey sediment with medium sand sized mica flakes", from 
ene tidal flats of Seal Beach Lagoon. Tables III through V 
contain detailed analysis of its mechanical properties, and 
Mineralogical and chemical analyses of the lagoonal mud from 
memples taken from the borrow pit. A year after the pene- 


tration tests were completed, Muga [1966], while measuring 
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TABLE Il. MECHANICAL PROPERTIES OF THE SEDIMENT IN THE TEST PIT 


Depth (in) 0-3 4-7 
Color (GSA No.) 5GY 6/1 5GY 6/1 
Odor HAS HS 
Bulk Wet Density (gm/cc) 1.452 lesog 
Vane Shear Strength (psi) 0.364 1.495 
Remolded Strength (psi) 0.364 1.736 
Sensitivity 1.0 0.7 
Water Content 82.9 66.3 
Specific Gravity of Solids 2.914 3,019 
Dry Density (gm/cc) 0.794 0.937 
Void Ratio | 2.676 2,226 
Porosity (%) 72.8 69.0 
Saturated Void Ratio 2.416 2.002 
Organic Carbon (%) 1.685 0.578 
Percent Sand 6.0 2.0 
Percent Silt 30.0 18.0 
Percent Clay 64.0 80.0 
Median Diameter (mm) 0.0020 0.0018 


Note: Samples were obtained with a two and one-half inch core liner from 
the site where the sediment originated. 


Seemples collected. by -R. WJ. ~Smith andjanalyzed at the 
Naval Civil Engineering Laboratory in Port Hueneme. 
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TABLE IV. MINERALOGICAL ANALYSIS OF THE TEST PIT SEDIMENT® 


COARSE SILT TO SAND PERCENTAGE 
Quartz joe 
K-Feldspar Tle 
Calcite 3.36 
Opal 1.47 
Muscovite (illite) bio. 


Total = 39.45 percent 


ACTIVE CLAYS 


Kaolinite 16,31 
Mixed Layer 22 27 


Total = 35,58 percent 


INACTIVE SILT TO. CLAY SIZES 


Goethite 5.83 
Halite 4.53 
Organic 4,4] 


Total = 14.77 percent 


Note: The mixed layer was composed of 67 percent montmorillonite and 33 
percent chlorite. The muscovite was chiefly terrigenous (clastic) and occurred 
in coarse silt~sized grains. The analyses were made by a combination of x-ray 
diffraction, chemical, and petrographic methods. Infrared absorption analyses 
and differential thermal analyses were used to obtain supplemental data. The 
values shown are recomputed to 100 percent and are reliable to about +15 
percent. 


Sihelan, J. A., S Reporte ion Mineralcoy ef Ocean Botcom 
Beaiments, p. 10 unpublished report to the Naval Civil 
menigineering Laboratory, September, 1967. 
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TABLE V. CHEMICAL ANALYSIS OF TEST PIT SEDIMENT” 


OXIDE PERCENTAGE 
$10; 48.34 
Aly 0, 18.26 
FeO, 7.88 
MgO 2.98 
CaO 1,92 
Ko™, 2.74 
Na,O 3.02 
HO 2.78 
Aor Ae : 9.67 
Total 97 59 
Surigeot 73.01 


2 


oe. 














TABLE VI. RANGE OF VARIATION OF TEST PIT SEDIMENT 


"The main feature of the soil is its highly variable properties, being 
both nonhomogeneous and nonisotrophic. Some of the variations are indicated 
below. 


RANGE OF VARIATION (42 TESTS) 


MINIMUM MAXIMUM AVERAGE 


Vane Shear Strength (psi) 0.0 1.050 0.530 
Remolded Strength (psi) 0.0 0.861 0.390 
Original Water Content (%) 43.5 726 63.5 
Void Ratio. 1.104 2.126 1.710 
Porosity S225 68.0 Giiee 
Liquid Limit 41.0 E850 5G, | 
Plastic Limit 27S 43.4 34,2 
Liquidity Index 70 255 136 
Percent Sand ed 25.5 10.2 
Percent Silt 30.8 54.4 40.5 
Percent Clay 35.8 66.0 46.4 


"In addition, soil samples were subjected to numerous triaxial, permeability 
and viscosity tests. From these tests it was determined that; (i) the yield 
strength of the '‘soil' is approximately 0.5 psi; (ii) permeability is practically 
nil; (iii) cohesion, c, is almost zero; (iv) angle of internal friction is almost 
zero and (v) viscosity is independent of the rates of shear or deformation 
ie., Tpm)." 


Pinca, J. Breakout. Wonces, pliay MoS. | Nawvad, Cava 
Baginecring Laboratory Technical Note N-863, Port 
memeneme, September, 1966: 
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meeak out forces in the test pit; conducted more tests fo 
merity the mechanical properties of the sediment; Table 
mieis a summaxry.of Muga's. tests. -Lb.i1sS ebvieus+that the 
fagoonal mud was not uniform but characterized by variability, 
mat a major cause of the variability was undoubtedly the 
mesult of sampling disturbances and continual reworking. 
Muga also stated that the mean viscosity was 300 square feet 
m_eeeeecond. This value represented a series of measurements 
Peaken with a commercial rotating-type viscometer. He did not 
Meeempt to correlate the viscosity with the vane shear 
Mereength because of the "marked variability in the soil 
Mmeoperties". Since the sediment contained approximately 

mee percent silt to sand sizes, there was certainly some 
mee Of internal friction caused by these large sizes and 
Bits Was probably the reason Muga found a high viscosity 
Peeeurement. It also appeared that there were more coarse 
meer sizes present a year after penetration tests were con- 
Gucted which indicated possible contamination with other 
Material. For this reason Muga's values were not taken as 
mEereaccly correct, but they were included to illustrate 
acta bliity of the sediment used for the tests. 

Pigure 12 Showsisthe mean shear strength profile of 
the sediment from 550 measurements. Although the mean 
merength profile is not representative of the average Shear 
Strength profile for the world's deep-ocean basins shown 
Seemed SOlid line, Figure 12, from Hoag [1970], it is 


remarkably close to the shear strength profiles of the edges 
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DEPTH (inches) 


SHEAR STRENGTH (psi) 
1 2 3 4 








mean 
one standard deviation from mean 





- average shear strength from 
114 cores from water depths 
greater than 6,000 
[Hoag, 1970] 





Figure 12. TEST PIT MEAN SHEAR STRENGTHS 
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of the Zigsbee Abyssal Plain as Shown by Bryant and others 
i967). 

The numbering system designating the location, or 
Beation, where each penetration test was. conducted is shown 
Sierigure 13. The shear strengths at Gach station were, taken 
meang a two inch by five inch vane. The vane was operated by 
meng £rom a punt floating on the surface. [It -was driven into 
See sediment in six inch increments and torque applied by 
hand with a standard torque wrench and the maximum torque, 
meaverted to shear strength, recorded. Not all of the vane 
Bests were conducted in the presence of an engineer and some 
mere done by different personnel. Generally the shear 
mepength profiles are consistently Similar, with some erratic 


memes, but there is some correlation between very erratic 


_eeiaes and the date of tests, which indicated operator 


“skill and care may have had ancwiniluence. The methodr or 


Obtaining shear strength values greatly influences the value 


Meeotained according to Terzaghi [1960] and it was difficult to 


~ sess 


get consistent readings from an unstable punt with no control 


meee Vane's rotational speed. 


Other methods of determining shear strength profiles 


m_emene mud pit were considered, but rejected for the follow- 


ang reasons. In 1965, the existing penetrometer probes 
related resistance encountered to shear strength by empirical 
tables or graphs, and the very low strength of the mud pit 
mes Out of the empirical range. The use of these devices 


would have required modifications of the instruments and 
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NOLE Small circle andueares the Locatxon of 
sub-lettered shear strength data points 
around a main station. 


Figure 13. LOCATION OF DROP STATIONS IN THE TEST PIT 
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extensive testing to obtain empirical relationships. Second, 
Ehe sediment in the mud pit probably had some slight but 
Mnknown angle of internal friction and ithe vane test mea= 
sures a cohesive shear strength which is influenced by the 
amegle of internal friction in a minor way. However, the 

mane test seemed the most practical test for shear strengths 
mucores and in-Situ testing. ‘Therefore, vane shear strengths 
mere used exclusively to keep the results comparable to a real 
Situation. 

The horizontal variability of the observed shear 
Seeenogths at six inch depths is shown in a-plan view of the 
Meeepit, Figure 14. The measurements taken on different 
Gates in the same location also show a Spread in values. 

Tie Variability Of thewshear strength in Ene mudypiaw 
Was probably not random for the following zeasons. 

om The huge: volume of Material requared tow) the 
m_eEemecessitated handling by large construction equipment 
meeen could have scooped up isolated chunks of material and 
metowed them to be placed in the pit essentially undisturbed, 
While most of the sediment was greatly disturbed during 
Meeavation, hauling and placement. The very high shear 
strengths have been attributed to the relatively small vane 
meoming in contact with one of these hard chunks. The chunk 
Meuld not necessarily affect penetration depth unless an 
Object was dropped very close to one because the sediment 
Seas incompressible and moved as a unit. Most of the objects 
meee Large compared to the chunks or vertical distance 


between anomalously high observed shear strengths. 
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NOTE: More than one value indicates measurements 
were taken on different dates. 


Figure 14. © SHEAR STRENGTHS AT SIX-INCH DEPTHS PLOTTED 
ON A PLAN VIEW OF THE TEST PIT (psi) 
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2. Although it was not noted in the analysis of 
the sediment, some organic material or foreign objects may 
have been present Since no special precautions were taken 
to Sieve the sediment while it was being placed. 

3. The extremely low walues of observed shear 
pesength may be attributed toe the method of handling the 
Bearment after a penetration trial. The holes made by the 
objects and other experiments with the sediment were smoothed 
eyer, filling depressions with a shovel and rake. Then shear 
merengths were taken-at—each station and penetration drops 
Menducted at a later time. After all the stations had-been 
Meeq, the procedure was repeated. Since other tests were 
meng conducted concurrently, the hole filling and smoothing 
Beene Continuous operation. Undoubtedly some very soft 
Material and water accumulated at, jche bottens of the holes 
ema this accounts for the anomalously low shear strengths at 
mepsiderable distances below the surface. 

tae low sensitivity o£, the sediment béetone Sampilaing 
ema the continual reworking it underwent while tests were in 
meegress, led to the conclusion that the sediment used for 
mee cests had a sensitivity of one and that the observed 
in-situ shear strength of the sediment was a good repre- 
SeueeetOn of its actual in-situ strength, if the problems of 
‘test Operator error and the anomalous points due to hard 
@hunks or soft spots could be resolved. 

Fortunately, enough data, 716 points, existed to make 


Meseatistical analysis and a mean strength profile of the 





| 
! 
. 
| 





moa pit. Table VIL includes a.matrix of el iethe ooserved 
Shear strengths. The anomalously high and low values were 
meaentified by eye. Since the profile should have been 
fairly smooth with all the reworking the sediment underwent 
meeing placement and tests, the hard and Soft Spots were 
SaoLtly selected by visual inspection. The values considered 
anomalous are underlined? 1m Tabile Wik. 

Tae mean and Sittandard deviation Of Sheam Strengths as 
a function of depth, both with and without the anomalous 
ees, are listed inTable VIII. It should be noted that 
eee anomalous values often cause the standard deviation to 
_meemenan double.-. The mean shear strength profile of the 
mae pit and one standard deviation calculated without the 
Smemalous points are shown in Figure 12. 

Me C€ach Station the Shontvous points were adjusted to 
meee a smooth profile that generally was within one standard 
Mevracion of the mean. The adjustment process was again 
Mm—mective, aS the precision of the original data did not 
mamtant a more sophisticated approach. The observed shear 
Btrengths for each station were plotted with the mean pro- 
file. Then the anomalous points were moved to fit a general 
Booth EEend OL points rani one standard deviation from the 
meen. The figures in Appendix B illustrate the adjustment at 
Gach station and the solid line is the actual shear strength 
profile used by the computer solution to predict penetration 
depth. 

The shear strength and standard deviation value at 


the air-sediment interface was taken to be zero, in that sea 
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TABLE VII. MEAN SHEAR STRENGTHS AND STANDARD DEVIATIONS OF 


TEST PIT, 
DEPTH MEAN SHEAR STRENGTHS ONE STANDARD DEVIATION 
(in) (psi) (psi) 

6" elles gebenacnes 23 38 
12" a ey hs 32 AA 
18" Om les 47 84 
24" 1.50 40) ieiies 63 1.21 
30" 1.80 4 eles 62 1.57 
36" 28 a Eas 70 1.94 
42" rer 2.87 43 60 1.87 
4g" ee Al 74 1.60 
54" ko 3.404 75 1.96 
60" 23s a5 7] 1.81 
66" 2.505 2.47 44 82 1.21 
72" 2.52 94 3.04 4 87 1.49 


Note: The first columns under the shear strengths and standard deviations are 
values computed with the anomalous points discarded. The second columns are 
values computed with all the data points, The subscripts in parenthesis are the 
number of data points. 
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Mater -was continually added to the tést pit. to offset 
Emaporation. losses and: the top of the pzt was’ continually 
Being ceworked- by raking. — In general, ajvirgin Undijsturbed 


meaiment Should have a finite strength at the surface. 


Go. RESULTS 

Table IX summarizes the observed penetration depths, the 
predicted depths, and tabulates the difference between ob- 
Served and predicted. A numerical analysis of the differences 
between observed and predicted depths was not made because, 
match few exceptions, the method developed accurately pre- 
@acted the penetration depth within the experimental uncer- 
SeaameLes Of depth measurement and shear strength accuracy. 
mie method overestimated penetration 17 times and under- 
meeamaced penetration 22 times. The largest differences were 
mmaerestimates of penetration. In all cases of differences 
Pee than four inches, the prediction had to rely on an 
me~istea Shear strength. . This indicated a problem with 
Measuring the shear strength at that station and implied the 
Seum@eemee Of a SOLL pocket in the sediment. The large dis- 
meemancies also occurred with impact velocities of five feet 
mer second. 

Maore X sets forth the penetratiqn data along with notes 
_werromn movies of the tests. The time-to-penetrate was 
Observed by counting the frames of the movie film from object 
impact until all downward motion had ceased. The rebound 
time was also observed by counting movie frames. The rebound 


Gistance was estimated from the pictures and is believed 
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Table IX. SUMMARY OF RESULTS 


OBJECT PENETRATION 
VELOCITY OBSERVED PREDICTED DIFFERENCE 
(ft/sec) (in) (in) (in) 
SPHERE (wt. = 56] Ibs) 
0.0 pe Loe7? ee 
0) hook, LA he = 13" 
10.0 14.0" 14.0" 0.0" 
15.0 242" 245" Ono 
20.9 20.0" 20.5" om Olen 
5,0 1620" jigs? - 4,5" 
10.0 16.03 12e le = 3.90 
fo .0 1730" 260" TAM 
20.0 19a 20.4" O90 
HORIZONTAL CYLINDER (wt. = 540 Ibs) 
0.0 508 6,0" Os 
3.0 Lives Lior - 4,0" 
10.0 13.0" 13.6" +0 38" 
15.0 17 70% 170" 0.0" 
20.0 1620" 1a.9" - 0.1" 
20.0 1422 14,7" 020" 
CYLINDER ON FLAT END (wt. = 540 Ibs) 
a0 14.0" S72 = 5.9" 
10.0 15250 12,8" 257" 
15.0 lo20" 14,2" - 0.8" 
20.0 4 £2 2 2 ae ame W.2" ee 
0.0 8.0" g.4" +0.4" 
0 4.2" 4,4" 44052" 
10.0 18.0" 16.0" - 2.0" 
13:0 14.0" Wy is = eo 
t35;,0 1s0* 14.6" ~ 3.4" 
20.0 liye Be 16.8" col AAO A 
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TABLE IX. (CONTINUED) SUMMARY OF RESULTS 


OBJECT PENETRATION 
MELOCITY OBSERVED PREDICTED DIFFERENCE 
(ft/sec) (in) (in) (in) 
CYLINDER ON FLAT END (wt. = 561 Ibs) 
220 b yy VeeSeRossas= 2.6%) Seay oO eeeee 
10.0 14,0" 14.0" 0.,.0% 
1o.0 te25? 20.4" Palen 
20.0 130" hee +2.6" 


CUBE (wt. = 702 lbs) 


0.0 a5" 3.7" ~0.8" 
5.0 we 4.4" hee 
10.0 8.5" 9.7" #429 
15.0 13.0" 12.8" -0.2" 
20.0 17.0" 14.9" aa 
RECTANGLE (parallelopiped, wt. = 1014 Ibs) 

0.0 3.0" 3.2" +0.2" 
5.0 5.0" 7.0" +2.0" 
10.0 9.2" 9.4" +0.2" 
15.0 14/8" 14.7" ie 
15.0 15.5" 19.0" +3.5" 
20.0 21.6" 20.1" =a" 
CONE (wt. = 1127 Ibs) 

0.0 55.0" 50.0" -5.0" 

5.0 53.5" 51.0" = 5 
10.0 49.5" 54.0" +4,5" 
Or 65.0" 64.0" - 1.0" 
20.0 62.5" 59.0" ~325" 
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BecuraclLe Eo one inch.) Penetration depth Trerers to the 
meyect's final resting position, dnG-did noe exelude che 
mebound distance. The percent of ultimate bearing capacity 
mes calculated as; the ultimate bearing Capacity Of the 
meaiment at the depth the object came to rest, divided by 
Seewectual force applied ‘to the sediment by the object at 
meeeinal resting position. The shear Strength adjusted 
eolumn indicates whether or not an anomalous shear strength 
Malue was adjusted which directly affected the calculations 
Seepenctration prediction. 

The time-to-penetrate measurement was done in an inexact 
Manner, but is remarkably constant when the range of 
metocities and different object shapes are considered. Its 
Betative constancy in relation to the other parameters may 
indicate thae ik 2s related CO a unique basve property oz 
fae sediment. Schmid [1966] observed the same phenomena 
Maen testing small penetrometers with different impact 

metocities. 

| ore cebound time and distance were very difficult to 
Estimate due to the short duration and fact that splatter 
Bsometimes partially obscured the object. Some of the rebound 
@e the object, and apparent elastic behavior of the sediment 
Bes Originally thought to be caused by a buoyancy effect. 

| However the cavity formed by penetration did not collapse 
Until rebound was over and the object was essentially at 
rest. Subjectively it appears that rebound is a function of 


Mempetastic property of the sediment. 
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The percent of ultimate bearing capacity that was reached 
Bethe different objects seems to relate (fo impact velocity. 
Maeall the penetration tests the objects were left in place 
mor at least 24 hours and in some cases up toa@ght days and 
een checked to see_if they had settled further into the 
meeement. No further settlement was observed in any case, 
Meme impact velocity of five feet—per—second the objects 
mais stantly penetrated to a level where 300 percent of the 
Setimate bearing capacity was available to support the object. 
mes implies that—a structure being emplaced by a free-fall 
mecnnigque would only require a velocity of five feet=-per— 
Becond to generate the safety factor of three recommended in 
See mechanics in a Situation where the ultimate bearing 
Capacity equation is used. 

ms observed in the movies, any air entrapped under .the 
mebject had negligible effect on either rebound or pene- 

Meaction depth. Rounded objects such as the sphere and 
m@orizontal cylinder did not behave differently from the flat 
meereo Objects that did entrap some air. 

Paaeecang and pile up of the Sediment was Slight. in ali 
Bcases, even with the rectangular parallelepiped impacting at 

P20 meee per second, the crater lip was less than-2-inches 
; 
high. The sediment moved laterally away from the objects and 


enot up. 
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The shear strength profile is the average of four 





sets of data from the location where the cones were 


No motion pictures were available for the 


dropped. 
cone tests. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A rational, phenomenological, theory has been applied Eo 
mre problems of free-fall penetration into a simulated deep-— 


mea sediment. The theory correlated well with full scale 


_ tests. A more, complete solution will probably only be 















available after more is understood about the mechanical 


' properties and dynamic behavior ‘of deep-sea sediment. 


A. CONCLUSIONS 
mre, method predicts the penetration of free—fatling 

Bp jects into weak, Saturated sediments for an air-to-sediment 
“case within the accuracy of the present state-of-the-art 
meni gues for measuring the sediment mechanical properties. 
It Should work equally well in a water-to-sediment case with 
‘Some modifications. 
| ae ainpact duration time, from the instant of contact 
mmtit the object comes to rest, was observed to be rela- 
tively eonstant and may be a unique property of the dynamic 
behavior of a Sediment type. The duration time was rela- 
Bively constant even though the penetration distance changed, 
mummeerject Shapes were different; and the kinetic energy of 
Se falling objects ranged from 0 to 5,000 foot-pounds. 
| The method may be used with larger, heavier objects but 

BCaution must be used in extending the method for use with 
smaller objects since viscous drag and surface adhesion 


| effects may be important for smaller masses. 
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At high impact velocities the effects. of added mass and 
form drag cannot be neglected) in icomputing penetration deprh. 
Mne accuracy of the constants and variables in the standard 
Bearing Capacity equation have the most control over pre- 


mieted penetration. depths. 


B. RECOMMENDATIONS 


















) 
> 


ime. method should: ‘be~experimentally verviied) by) full 
mae ECSLS. at Sea or by controlled tests inia Large: tank 
ashore uSing a uniform sediment whose mechanical properties 
are mHOWn. »- EXpPeriments, ef the scale. needed would, requaxecva 
Major effort but could include investigations into other 
Munknown areas such as the free-fall attitude and stability 
ewer ious shapes in water, the verification of terminal 
Welocities of free-falling unstabilized objects in a. water 
Motumh and the study of breakout forces and Salvage methods 
mmmeed to retrieve the objects. 
Whenever an object is recovered from the deep-ocean floor 
@very effort should be made to obtain photographs of the 
mMbject. The sediment around the object should be extensively 
Bampled and completely analyzed to allow verification of this 


meemod Of predicting penetration depth. 
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APPENDIX A 


FORCES ON OBJECTS PENETRATING A TYPICAL SEDIMENT 


mppendix A includes a set Of ‘cunves@whten. qualia tatively 









Seeerustrate the relative influence of each force in the 
spenetration equation. They were obtained by using the 
meeument shear strength profile typical of the world's 
deep-ocean basins as determined by Hoag [1970] and shown 
meeeagure 13. The curves are normalized for depen Gf péne= 
tration on the abcissa and the maximum force encountered on 
the Becdinate. Time 1S.independent of the calculations and 
Bre Pepe Of the curve ‘near the origin jas largely a function 


the iteration interval. 
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IMPACT VELOCITY EQUALS 


SEDIMENT SHEAR STRENGTHS IN HALF FOOT INTERVALS FROM SURFACE ARE 
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IMPACT VELOCITY EQUALS 


SEDIMENT SHEAR STRENGTHS IN HALE FOOT INTERVALS FROM SURFACE ARE 
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IMPACT VELOCITY EQUALS 10.0 


SEDIMENT SHEAR STRENGTHS IN HALF FOOT INTERVALS FROM SURFACE ARE 
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WORK DONE IN SECTION — FOOT POUNDS 


CONDITIONS AT THE END OF SECTION 
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